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Abstract 

We extract the form factors relevant for semileptonic decays of D and B mesons from 
a relativistic computation on a fine lattice in the quenched approximation. The lattice 
spacing is a = 0.04 fm (corresponding to = 4.97 GeV), which allows us to run very 
close to the physical B meson mass, and to reduce the systematic errors associated with 
the extrapolation in terms of a heavy quark expansion. For decays of D and Df. mesons, 
our results for the physical form factors at = are as follows: f^^^{^) = 0.74(6)(4), 
/+^^(0) = 0.78(5)(4) and /+^^^(0) = 0.68(4)(3). Similarly, for B and Bs we find: 
/f^'^(O) = 0.27(7)(5), /f^-^(O) = 0.32(6)(6) and /+=^-^(0) = 0.23(5)(4). We compare 
our results with other quenched and unquenched lattice calculations, as well as with light- 
cone sum rule predictions, finding good agreement. 



PACS numbers: ll.15.Ha, 12.38. Gc, 13.20.Fc, 13.20.He 



1 Introduction 



Heavy meson decays are the main source of precision information on quark flavor mixing pa- 
rameters in the Standard ModeL The over-determination of the sides and the angles of the 
CKM unitarity triangle is the aim of an extensive experimental study: It addresses the question 
whether there is New Physics in flavor- changing processes and where it manifests itself. One of 
the sides of the unitarity triangle is given by the ratio Vet is known to approximately 

2% accuracy from b ciui transitions [1], [2] whereas the present error on Vub is much larger 
and there is also some tension between the determinations from inclusive and exclusive decay 
channels. Reduction of this error requires more experimental statistics but — even more so — an 
improvement of the theoretical prediction of the semileptonic spectra and decay widths. 

This is the prime motivation for the study of semileptonic form factors of decays of a heavy 
meson H = B,D into a light pseudoscalar meson P = 7r,K, which are usually defined as 

(PipW^lHipn)) = + (p'h + p' - '"^~'"^ ^) Mq') ■ (1) 

Here = q2l^(li is the vector current in which qi (^2) denotes a light (heavy) quark field; p 
{Ph) is the momentum of the light (heavy) meson with mass mp {mn), and q := pn — p is the 
four-momentum transfer. The fo{q^) and f+{q^) form factors are dimensionless, real functions 
of q^ (in the physical region), which encode the strong interaction effects. They are subject to 
the kinematic constraint /+(0) = /o(0). 

In the approximation of massless leptons (which is highly accurate for £ = e or £ = /i) , the 
differential decay rate for the H Pii^i process involves f+{q^) only: 



dq^ 1927r3m^ 



(m^ + ml-q^Y - Amlrnl] | /+ (g') | ' . (2) 



Another motivation for our study is that fo{q^) and f+{q^) enter as ingredients in the analysis 
of nonleptonic two-body decays like B —>■ irn and B — > irK in the framework of QCD factor- 
ization [3l H], with the objective to extract CP- violating effects and in particular the angle a 
of the CKM triangle. One issue that is especially important in this respect is the question of 
flavor SU(3) violation in the form factors of the decay B ^ ir vs. the rare decay B K. 

High-statistics unquenched lattice calculations of D-meson (and also S-meson) decay form 
factors in the kinematic region where the outgoing light hadron carries little energy (small 
recoil region) have been performed recently [5], [6l [7] and attracted a lot of attention. Direct 
simulations at large recoil, q^ <^ m^, with light hadrons carrying large momentum of order 
2 GeV, prove to be difficult and require a very fine lattice which is so far not accessible in 
calculations with dynamical fermions. This problem is aggravated by the challenge to consider 
heavy quarks which either calls for using effective heavy quark theory methods or, again, a 
very fine lattice. In practice, one is forced to rely on extrapolations from larger momentum 
transfer q^ and/or smaller heavy quark masses. Several extrapolation procedures have been 
suggested [H [H [TOl [11], [12] that incorporate constraints from unitarity and the scaling laws in 
the heavy quark limit. Alternatively, 5-meson form factors in the region of large recoil have 
been estimated using light-cone sum rules [121 HI] (for recent updates see [151 [ISl [13 [IB])- 

In this paper we report on a quenched calculation of semileptonic H —>■ Piue form factors 
with lattice spacing a ~ 0.04 fm using nonperturbatively 0(a) improved Wilson fermions and 
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0{a) improved currents. On such a fine lattice a relativistic treatment of the c quark is justified 
and also the extrapolation to the physical b quark mass becomes much more reliable compared 
to similar calculations on coarser lattices. In addition, we can explore possible subtleties in 
approaching the continuum limit in form factor calculations: In our previous work [19] we did 
find indications for a substantial discretization error in the decay constants /^^ etc.; similar 
conclusions have also been reached in Ref. [20]. This is particularly relevant in view of the 
claims of evidence for New Physics from comparison with recent dynamical simulations — see, 
e.g. Ref. [21] for a discussion. 

On physical grounds, one may expect a nontrivial continuum limit because form factors at 
large momentum transfer are determined by the overlap of very specific kinematic regions in 
hadron wave functions (either soft end-point, or small transverse separation). The common 
wisdom that hadron structure is very "smooth" — and that numerical simulations on a coarse 
lattice could thus be sufficient to capture the continuum physics — may not work in this par- 
ticular case. This can be especially important for SU(3) flavor-violating effects, which are of 
major interest for the phenomenology. Inclusion of dynamical fermions and the approach to the 
chiral limit are certainly also relevant problems, but not all issues can be addressed presently 
within one calculation. 

This work should be viewed as a direct extension of the investigation of the APE collaboration 
in Ref. [22], who performed a quenched calculation with the same nonperturbatively 0{a) 
improved action and currents. Also their data analysis is similar. However, they use coarser 
lattices with a ~ 0.07 fm [f] = 6.2). On the other hand, the spatial volume of their lattices is 
very close to ours (L ~ 1.7 fm). So the main difference lies in the lattice spacing, and a direct 
comparison of the results is possible yielding information on the size of lattice artefacts, while 
there is no need for us to perform simulations on a coarser lattice ourselves. 

The presentation is organized as follows. Our strategy is discussed in detail in Sec. [21 It 
allows us to run fully relativistic simulations for values of rriH up to the vicinity oims'- This is 
achieved by using a lattice characterized by a very fine spacing a. The extraction of physical 
quantities from our data and the final results with the associated error budget are presented 
in Sec. El The final Sec. H] contains a summary and some concluding remarks. Some technical 
details and intermediate results of our calculation are shown in the Appendix. Preliminary 
results of this study have been presented in Refs. [231 12^ - 

2 Simulation details 

The lattice study of heavy hadrons is an issue that involves some delicate technical aspects: 
The origin of the problem stems from the fact that, typically, the lattice cutoff is (much) smaller 
than the mass of the B meson. 

Common strategies to solve this problem are based on heavy quark effective theory (HQET), 
i.e. expanding the relativistic theory in terms of where mq is the mass of the heavy quark. 
One can simulate in the static limit [25] or keep correction terms in the action to simulate 
at finite mq (non- relativistic QCD or NRQCD) [26]. These approaches have been employed 
effectively for studying B physics (see, for example, Refs. [271 [281 l29l [30] ). 

However, for smaller quark masses like the c quark in D mesons a large number of terms in 
the expansion must be included, making the simulations less attractive. The Fermilab group 
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developed an approach which interpolates between the heavy- and light-quark regimes [3T] . 
The coefficients accompanying each term in the action are functions of the quark mass and 
in practice, normally, the lowest-level action is used. This corresponds to using the 0(a) 
improved relativistic action with a re-interpretation of the results. Except for HQET [32l [33] , 
the associated renormalization constants for these approaches are only known perturbatively. 

We reduce the uncertainties related to the extrapolation to the physical heavy meson mass 
by using lattices with a small spacing in conjunction with a non-perturbatively improved 0{a) 
relativistic quark action. This theoretically clean approach enables one to get sufficiently close 
to the mass of the physical B meson, so that the heavy-quark extrapolation is short-ranged. 
In addition, in the region of the D meson mass, the discretization errors are reduced to around 
1%, see Section ISTTl 



L^xT 


403 X 80 




6.6 


lattice spacing a 


0.04 fm 


physical hypervolume 


1.6^ X 3.2 fm^ 




4.97 GeV 


^ of configurations 


114 


^critical 


0.135472(11) 


l^heavy 


0.13, 0.129, 0.121, 0.115 


I^Hght 


0.13519, 0.13498, 0.13472 


mp 


526 MeV, 690 MeV, 856 MeV 


csw 


1.467 


Zv 


0.8118 


hv 


1.356 


Cv 


-0.0874 



Table 1: Parameters of the lattice calculation (see the text for the definition of the various 
quantities) . 

Table [1] summarizes basic technical information about our study. We use the standard Wilson 
gauge action to generate quenched configurations with the coupling parameter j3 = 6.6. For this 
parameter choice, the lattice spacing in physical units determined from Ref. [51] using Sommer's 
parameter tq = 0.5 fm is a = 0.04 fm. Our calculation is based on the 0{a) improved clover 
formulation for the quark fields [35], with the nonperturbative value of the clover coefficient csw 
taken from Ref. [36]. We use 0{a) improved definitions of the vector currents in the form [37] 

(^27m9i + iacvd^q2a^^qi) (3) 

with = I [7/i,7i/]. The renormalization factor Zy, the improvement coefficient by as well 
as Cv are known nonperturbatively [381 EHl SO]- AH statistical errors are evaluated through 
a bootstrap procedure with 500 bootstrap samples. We consider three hopping parameters 
corresponding to "light" quarks, hiught (the corresponding masses of the light pseudoscalar 
meson states mp are also given in Table [1]), and four hopping parameters, i^heavy, corresponding 
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Figure 1: Diagram representing the semileptonic decay of a heavy-light pseudoscalar meson 
to a hght pseudoscalar meson (left panel). A schematic representation of the corresponding 
three-point correlator calculated on the lattice is also shown (right panel). 

to "heavy" quarks; in particular, n = 0.13498 and k, = 0.129 are found to correspond to quark 
masses close to the physical strange and charm quark mass, respectively. 

The extraction of the matrix element appearing in Eq. ([T]) from the lattice can be done 
by considering the large time behavior of three-point correlation functions {0,tx,ty) for a 
pseudoscalar light meson sink at time t = 0, a vector current at time tx, and a pseudoscalar 
heavy-light meson source at time ty = T/2 (see Fig. [1]): 

Cf)(0,t.,t,) = 5^e-P^--V^"^^(i7^(y,t,)V;(f,t,.)^^(0)) . (4) 

Here, and are Jacobi-smeared operators of the form q^l^qs and qi'j5qs, respectively; qt 
denotes the heavy quark, qi is the decay-product quark, while is the "spectator" quark. 

For sufficiently large time separations (i.e. ^ t^,- ^ T/2 or T/2 ^ t^; <^ T), cjp {0,tx,ty) 
behaves as: 

f S5e-^*^e-^-(*-*^) {H{ph)\V,\P{p)) for < T/2 

Cf(0,t.,g^<^ , (5) 

[ ±lll^e-^(^-'^^e-^"('^-'yUH{pH)\V,\P{p)) for > T/2 

with Zfj = \{0\H^\H{ph))\ and = |(0|P"^|P(p)) |, while E (Eh) denotes the energy of the 
light (heavy) meson. To extract the matrix elements we divide the three-point functions by the 
prefactors, which are extracted from fits to smeared-smeared two-point functions. The matrix 
element is then obtained by fitting this result to a constant, in an appropriate time range where 
a clear plateau forms (for example, for 12 < < 28). 

We consider three-point functions associated with different combinations of the momenta p 
and ph, which are listed in Table [21 In particular, we focus our attention onto three-momenta 
of modulus and 1 [in units of 27r/(aL)], since they yield the most precise signal, restricting 
ourselves to the cases where p and pn lie in the same direction. Thus we measure directly 
5 different values for the form factors, for every Kught and Kheavy combination. The full form 
factors can then be constructed from the data points obtained this way, by making an ansatz 
for the functional form of fo{q^) and f+{q^). 

In the present work, we fit our data with the parametrization proposed by Becirevic and 
Kaidalov [S]: 

1 — 1 — (7^ 1 — ao^ 
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Ph 


P 




(0,0,0) 


(1,0,0) 


(-1,0,0) 


(1,0,0) 


(-1,0,0) 


(2,0,0) 


(0,0,0) 


(0,0,0) 


(0,0,0) 


(1,0,0) 


(0,0,0) 


(1,0,0) 


(1,0,0) 


(1,0,0) 


(0,0,0) 



Table 2: Momentum combinations considered in the analysis of the three-point functions, in 
units of l-R j[aV). 

where q := q/mii*, ttih* being the mass of the lightest heavy-light vector meson. 

The parametrization for the form factors given in Eq. ([6]) accounts for the basic properties 
that come from the heavy-quark scaling laws in the limits of large and small recoil and also 
satisfies the proportionality relation derived in Ref. [H]. It is also consistent with the trivial 
requirement that the l.h.s. of Eq. ([T]) be finite for vanishing momentum transfer, which implies 
/o(0) = /+(0). The results that we obtained for the three parameters entering Eq. ([6]) from a 
simultaneous fit to /o and /+ are presented in the Appendix. 

Some alternative ansatze for the functional form of f+{q^) were proposed in Refs. [IS, [9l [10] 
and are discussed in Ref. [12]: They yield results essentially compatible with each other and 
with the Becirevic-Kaidalov parametrization Eq. ([6]). More recently, Bourrely, Caprini and 
Lellouch [12] discussed the representation of f+iq"^) as a (truncated) power series in terms of 
an auxiliary variable z. A similar parametrization has also been recently used by the Fermilab 
Lattice and MILC collaborations, see Refs. [l3l HI] for a discussion. 

3 Extraction of physical results 

In order to extract physical results from our simulations, we follow a method analogous to 
Ref. [22]. We first perform a chiral extrapolation in the light quark masses. For a given 
quantity $ [one of the BK parameters appearing in Eq. (I6l)], the extrapolation relevant for 
decays to a pion is performed as follows: We fit the results obtained at different values of the 
mass of the pseudoscalar state linearly in m^, 

$ = Co + ci ■ , (7) 

and extrapolate to mp = m^, where m,r is the mass of the physical pion. Examples of the 
extrapolations are shown in Figure [2] for the case of $ = /+(0), a and j3 at Hheavy = 0.115. On 
the other hand, for decays to a kaon, we hold the hopping parameter of one of the two final 
quarks fixed to k, = 0.13498, which, for our configurations, corresponds to the physical strange 
quark at a high level of precision [IS], and perform a short-ranged extrapolation of the curve 
obtained from the linear fit in mp to the square of the mass of the physical K meson. 

Then we perform the interpolation to the physical c quark mass in terms of a heavy-quark 
expansion for the D (or Dg) meson decays, or the extrapolation to the physical b quark mass 
for the B (or Bg) meson. For our data, the extrapolation of the heavy quark mass to the 
physical b mass is short-ranged: for the heaviest Hheavy = 0.115, it turns out that the inverse 
of the pseudoscalar meson mass (with the light quark mass already chirally extrapolated to its 
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Figure 2: Extrapolation of the Becirevic-Kaidalov parameters to the chiral hmit, for decays to 
a pion, at a fixed value i^heavy = 0.115. The parameters obtained for decay product = f^spectator 
are extrapolated linearly in mp. The extrapolated values are shown as the full black dots. 



Decay 


fit 


lo 


h 


h 


xVd.o.f. 


B,D ^71 


linear 
quadratic 


^•-L-l.o 
S 1+2-9 


-4.11^1 

_q Q+6.9 


K q+8-5 


0.1377/2 
0.021/1 


B,D^ K 


linear 
quadratic 


^•^-0.9 


'-'■^-1.9 
-| 9 9+6.3 
-"-^■^-8.2 


7 7+7.3 
' • ' -5.4 


0.3247/2 
0.03813/1 




linear 
quadratic 




-2 

-11 n+''-'^ 


q 7+5-4 
' -4.6 


0.4025/2 
0.001888/1 



Table 3: The coefficients obtained from the fits to /+(0) in powers of according to 
Eq. ([8]) for different decays. 

physical value) is about = 0.243 GeV~^, to be compared with m~^^ = 0.189 GeV~^ for the 
physical B meson. The extrapolation can be performed by taking advantage of the fact that, in 
the infinitely heavy quark limit, the Becirevic-Kaidalov parameters appearing in Eq. ([6]) enjoy 
certain scaling relations: cbk\/it^h^ {P — ^)'^h and {l — a)mH are expected to become constant 
in the tuh — ^ oo limit. For finite one can parametrize the scaling deviations in powers of 
m-jl: 

ip = Iq + Ii- + I2 ■ + . . . (8) 

where ip G {cbk^/tt^h^ (/^ ~ 1) ■ '"^Hi (1 — «) ■ rnu}. Note that, since /+(0) = cbk ■ (1 — one 
can also use = /+(0) ■ m'^^ — which was, in fact, our choice. 

The extrapolation of m^^/+(0) is presented in Figure [31 The figure clearly shows the ad- 
vantage of simulating on a fine lattice, which allows us to probe a mass range very close to 
the physical B meson mass. We compare the results obtained from an extrapolation to the 
inverse of the physical B meson mass using either a first- or a second-order polynomial in mj^^ 
for the fit function, finding consistency (within error bars), for all decays. The corresponding 
fit results are listed in Table [31 In the following we refer to this first method as the "coefficient 
extrapolation" method. 

An alternative method to extract the physical form factors from the lattice data was pro- 
posed by the UKQCD collaboration [IS]. It consists of performing the chiral and heavy quark 
extrapolations at fixed v - p = {m^^ + m?p — q^)/ (2m//) , where v is the four- velocity of the heavy 
meson and p is the four-momentum of the light meson. The following steps are performed: 

1. fit of the form factors measured from the lattice simulations to the parametrization in 



6 



B, Dton decays 



y 4- 
> 

as- 



S 1 




B , D to K decays 



s s 

n ' r 



> 

o 




1- 



0.1 



0.2 0.3 0.4 0.5 



0.6 0.7 



a 



0.1 



0.2 0.3 0.4 0.5 

[GeV"^] 



0.6 0.7 



Figure 3: Left panel: The green (red) squares denote the interpolated (extrapolated) form 
factor mJ'V+(0) to the physical D (B) meson, for a decay to a pion, using a linear fit in l/rriH 
(solid line). A quadratic fit to the data is also shown (dashed line). Right panel: the results 
for the case of decay of a Dg (Bg) meson into a kaon. 



Eq. dSD; 

2. interpolation of the form factors at given values of v ■ p within the range of simulated 
data; 

3. chiral extrapolation of the points thus obtained, via a linear extrapolation in mp to either 
ml or mj^ (as described above); 

4. linear or quadratic extrapolation in mjj^ to the inverse of the physical heavy meson mass 
for the quantities: 



70 
2/3o 



as{mH)_ 



f u 



(9) 



which enjoy scaling relations at fixed v-p [46J|47]. Here, Po is the first /3-function coefficient, 
while 7o = —4 denotes the leading-order coefficient of the anomalous dimension for the 
vector current in HQET. It yields a (subleading) logarithmic dependence on rriH — see 
also Refs. |22l US] for further details; 



5. final fit of the points thus obtained to the parameterization in Eq. 

For comparison, we also calculate the physical form factors using this alternative approach, 
finding consistent results. This is illustrated in Table H] which summarizes the results for /+(0) 
from both methods. 

For our final results we take those obtained from the coefficient extrapolation method. We 
found this method to be superior in our case as the UKQCD method suffered from the fact that 
there was only a small region of overlap in the ranges of v ■ p for the form factors at different 
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Klight and Kheavy In addition, since the data can be fitted with both a hnear and quadratic 
function in mjj^, we use the hnear fits for the central values and statistical errors and use the 
differences in the results from the linear and quadratic fit to estimate the systematic errors, as 
discussed in the next section. Our results for the form factors at finite are shown in Figs. H] 
and [51 





Coefficient 


extrapolation 


UKQCD method 


Decay 


linear in m^^ 


quadratic in mjj^ 


linear in m^^ 


quadratic in m^^ 


D ^TT 


0.74+[! 


0.731^ 




0.69+^ 


D ^ K 


0.781^ 


0.77+^ 


0.75+^ 


0.75+^ 


Ds^K 


0.68l^ 


0.671^ 


0.68+^ 


0.67+^ 


B 


0.271^ 


0.30^^^ 


0.29+^^ 


0.3111^ 


B ^ K 


0.321^ 


0.35^^ 


0.35t^i 




Bs ^ K 


0.231^ 


0.26+^ 


0.23+^ 


0.271^ 



Table 4: Final results for the physical values of the /+(0) form factor, for different decays, 
with statistical errors only. We compare the results obtained from the coefficient extrapolation 
and UKQCD methods as well as different truncations of the heavy quark expansion when 
extrapolating or interpolating in m~^. 



3.1 Systematic uncertainties 

Systematic uncertainties affecting our lattice calculation include: the quenched approximation, 
the method to set the quark masses, the chiral extrapolation for the light quarks, discretization 
effects, the extrapolation (interpolation) of the heavy quark to the physical h (c) mass, finite 
volume effects, uncertainties in the renormalization coefficients, and effects related to the model 
dependence for /o,+ (g^). Let us now consider each source of error in turn. 

Quenched approximation: the size of the error this approximation introduces is not 
known. However, one can take as an estimate the variation in the results if different quantities 
are used to set the scale. In the quenched approximation different determinations of the lattice 
spacing vary by approximately 10% [H]. By repeating the full analysis, we find that varying 
the lattice spacing by 10% induces an uncertainty of approximately 2% for the D n decay, 
and of approximately 12% for i? — *■ vr. 

Setting the quark masses: we use the k values corresponding to the light (u/d) and 
strange quarks determined in Ref. [19]: ki = 0.135456(10) and Hg = 0.134981(9). The uncer- 
tainty in these determinations leads to a very small uncertainty in the form factors. For the 
c and b quarks we do not quote the corresponding k values. We interpolate (or extrapolate) 
our results directly to the physical masses of the pseudoscalar heavy-light states. The result- 
ing uncertainty is determined by the statistical errors of the masses used for the interpolation 
(or extrapolation). The latter are found to contribute only a negligible amount to the overall 
systematic uncertainty. 

Chiral extrapolation: the method we used to perform the chiral extrapolation of our 
simulation results is discussed above. Note that the use of a large lattice practically constrains 
us to use only a few and relatively large values for the light quark mass (so that the masses 
of our lightest pseudoscalar mesons are far from the physical pion mass). However, as the 
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Figure 4: Physical form factors for D and Dg decays as a function of from this work and 
other quenched and dynamical studies. The solid black lines are the form factors obtained from 
the coefficient extrapolation method where Eq. ([8]) has been truncated at 0(mj^^), while the 
dashed black lines indicate the error on the form factors. The range oi v ■ p values achieved in 
our simulations approximately corresponds to —1.5 GeV^ ^ Q'^ ^ 2 GeV^. The dashed red lines 
are the results for the coefficient extrapolation method from Ref. [22]. The open red squares 
and circles are their results obtained using the UKQCD method. 

examples in Figure [2] show, the dependence of our results on the light quark mass is rather 
mild. So the size of the uncertainty arising from the chiral extrapolation though difficult to 
estimate is unlikely to be large. 

Discretization effects: as it was already remarked above, the leading discretization effects 
in our calculation are reduced to O(a^); given that our lattice is very fine (a = 0.04 fm), the 
associated systematic error can be estimated to be of the order of 1% (10%) for the decays of 
charmed (beautiful) mesons |19j . 

Extrapolation/interpolation of the heavy quark: our data can be fitted to both a 
linear and quadratic function in m'j^ with a reasonable x^. We use the results for the linear 
fit for our final results and the difference between the linear and quadratic fit as an indication 
of the systematic uncertainty. This leads to approximately a 1% uncertainty for D decays and 



9 





5 10 15 20 25 5 10 15 20 

[GeV^] [GeV^] 



Figure 5: Same as in Fig. HJ but for B and Bg decays; in this case, the v ■ p values of our 
simulations are in the range 14 GeV^ ^ ^ 23 GeV^. For B ^ n, the dashed and solid 
magenta lines in the range = — 14 GeV^ indicate the prediction from light-cone sum 
rules [I51[I6]. 



8% uncertainty for B decays. 

Finite volume effects: for our calculation, finite-volume effects are not expected to be 
severe; in particular, the correlation length associated with the lightest pseudoscalar state that 
we simulated (for nught = 0.13519) corresponds to approximately 9 lattice spacings, which is 
more than four times shorter than the spatial extent of our lattice. Systematic infrared effects 
can thus be quantified around 1-2%. This is comparable with the estimate of Ref. [H], in 
which, using chiral perturbation theory |^ [50] , the finite volume effects for their calculation 
with 2-1-1 flavors of staggered quarks and values of mpL between 4 and 6 are estimated to be 
less than 1%. 

Renormalization coefRcients: the uncertainty associated with the Zy coefficient, as de- 
termined in Ref. [38] for the quenched case, is about 0.5%. The same article also quotes a 
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1% uncertainty for by, which induces an error about 1% for decays of D mesons and about 
3% for B mesons. Concerning cy, a look at the results displayed in Fig. 2 of Ref. [39] would 
suggest that the relative error in the region of interest {g^ ^ 0.91) may be quite large, around 
30%; however, it should be noted that cy itself is a relatively small number, of the order of 
9%, and the impact of the uncertainty on cy on our results is about 1% (2%) for decays of D 
(respectively: B) mesons. 

Model dependence: finally, the systematic effect related to the ansatz to parametrize the 
form factors was estimated in Ref. [12], through a comparison of different functional forms that 
satisfy analogous physical requirements. For the B ^ ir decay, it turns out to be of the order 
of 2%. 

Combining the systematic errors in quadrature we arrive at an overall error of 5% for D 
decays and about 18% for B decays. 

3.2 Comparison with previous results 

Our results can be compared to other lattice calculations of these quantities and also with 
results of light-cone sum rules (LCSR) [13], [Hj. Table [5] summarizes the comparison for /+(0), 
while for finite the form factors from other studies are displayed in Figs. H] and [5l In the 
following we discuss in detail the comparison with these works, highlighting the advantages and 
limitations of the different approaches, as well as the possible sources of discrepancies. 

Our results can be closely compared with those obtained by the APE collaboration in 
Ref. [22], reporting a calculation very similar to ours. They worked in the quenched approx- 
imation, using the same non-pertubatively 0(a) improved action and currents and a similar 
analysis; on the other hand, their simulations were performed on a coarser lattice, with jS = 6.2, 
yielding a lattice spacing a = 0.07 fm, or a^^ ~ 2.7 GeV. The table and figures show that their 
values for the form factors lie around 3a {D tt) and 2.5a [D — > K) below our results, in terms 
of the statistical errors, in the region of = 0. If we adjust the APE results to be consistent 
with setting the lattice spacing using tq instead of the mass of the K* (used in Ref. [22]), the 
discrepancy reduces slightly, down to roughly 2.5cr {D tt) and 2a {D K). Assuming that 
0{a?) errors are the dominant source of the discrepancy, the difference in the results of the two 
studies is consistent with an upper limit on the discretization errors of approximately 0.08, or 
slightly above la in our results for /+^''(0) and 0.23 or 3 - 4(T in the APE results. 

For B decays we are not able to make such a close comparison, because the study in Ref. [22] 
extrapolates to the B meson from results in the region of 1.7 — 2.6 GeV for the heavy-light 
pseudoscalar meson mass. Although one would expect larger discretization effects for the 
B decay form factors, we find close agreement between our values and those from the APE 
collaboration. However, we should point out that any potential discrepancy may be masked by 
the long-ranged extrapolation in the heavy quark mass. 

Several unquenched calculations have been performed recently, which are based on the MILC 
Nf = 2 + 1 dynamical rooted staggered fermions configurations [51]. Results are available from 
joint works from the Fermilab, MILC and HPQCD collaborations for D decays [5], and from 
Fermilab and MILC [52l 144] and (separately) HPQCD [7] for B decays. These results were 
obtained using the MILC "coarse" lattices with a = 0.12 fm for D decays and including a finer 
lattice with a = 0.09 fm for the B decays. While these lattices are much coarser than those 
used in both our and the APE study a detailed analysis of the chiral extrapolation was possible 
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through the use of 5 hght quark masses for the 0.12 fm lattice (only two values were used for 
a = 0.09 fm). 

The Fermilab, MILC and HPQCD joint work for D ^ vr and D ^ K used an improved 
staggered quark action ("Asqtad") [51] for the light quarks and the Fermilab action for the 
heavy quark. To the order implemented in the study, the Fermilab action corresponds to a re- 
interpretation of the clover action. This approach can be used to simulate directly at the charm 
and bottom quark mass at the expense of more complicated discretisation effects. Discretisation 
errors arising from the final state energy (5%) and the heavy quark (7%) are estimated to lead 
to the largest systematic uncertainties in the calculation (compared to the 3% error from the 
chiral extrapolation). Given the coarseness of the lattice used, repeating the analysis on a 
much finer lattice would enable the estimates of the the discretisation errors to be confirmed. 
Overall, the results are consistent with ours, which suggests that the systematic effects due to 
the quenched approximation are not the dominant source of error. 

For the decay B ^ ir, Fermilab and MILC used the same quark actions as for the study 
of D decays. Using the 5 light quark masses at a = 0.12 fm and 2 light quark masses at 
a = 0.09 fm they performed a joint continuum and chiral extrapolation which removed some 
of the discretisation effects. They estimated that a 3% discretisation error arising from the 
heavy quark remains after the extrapolation. The results at finite are compared with ours 
in Figure O with statistical and chiral extrapolation errors only (which cannot be separated). 
A value for /+(0) is not given in Ref. [H] which focuses on extracting \Vub\ at finite using 
the parameterisation of Bourrely, Caprini and Lellouch [12]. However, an earlier analysis on 
the 0.12 fm lattices only was reported in Ref. [52]. Their result for /+(0) is given in Table [5l 

HPQCD performed the calculations for the i? — *■ tt decay on the MILC configurations using 
Asqtad light quarks and NRQCD for the b quark. Use of the latter enables direct simulations 
at the b quark mass. However, as NRQCD is an effective theory the continuum limit cannot 
be taken and scaling in the lattice results must be demonstrated at finite a. Results from the 
coarse lattice are shown in Figure O with statistical and chiral extrapolation errors only and for 
/+(0) in Table [5l A limited comparison of results on the finer lattice for one light quark mass 
did not indicate that the discretisation errors are large. The systematic errors are dominated 
by the estimated 9% uncertainty in the renormalisation factors which are calculated to 2 loops 
in perturbation theory. 

The Fermilab-MILC and HPQCD results are consistent with each other to within 2a and 
are also consistent with our results and those of the APE collaboration. As for the studies of 
D decays this suggests that quenching is not the dominant systematic error in the calculation 
of B ^ 71 decay. Similarly, unquenched results on finer lattices are needed to investigate the 
discretisation effects. Finally, note that in order not to overload Figures H] and [5], we do not 
show the (older) quenched results of the Fermilab group [53] . For i? — *• tt decays these results 
are within the range of the other existing calculations, whereas for D-decays the form factors 
come out 10 — 20% larger compared to most other calculations and also the new unquenched 
results obtained with similar methods. 

A different type of comparison can be made with the estimates obtained in the framework of 
LCSR. This analytical approach is, to some extent, complementary to lattice calculations, since 
it allows one to calculate the form factors directly at large recoil, albeit with some assumptions. 
Figure O compares our extrapolation of the f^^^{q^) form factor in the region < 12 GeV^ 
with the direct LCSR calculation [T5l [T6] . Their predictions are compatible with our results. 
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Similar consistency is found between lattice and LCSR calculations of /+(0), as seen in Tabled 
for both B and D decays. Note that the uncertainty quoted for /+(0) for B decays is smaller 
than that for D meson decays, and comparable with the precision of the lattice results. However, 
while LCSR provides a systematic approach for calculating these quantities it is by definition 
approximate and the errors cannot be reduced below 10 — 15%, unlike the lattice approach, 
which is systematically improveable. 



Decay 


This work 


Other results 


Source 


Method 


D 


0.74(6)(4) 


0.64(3)(6) 
0.57(6)(1) 
0.65(11) 
0.63(11) 


Fermilab-MILC-HPQCD 
APE |22] 
Khodjamirian et al. [M] 
Ball [Sg 


Nf = 2 + 1 LQCD 
iV/ = LQCD 
LCSR 
LCSR 


D ^ K 


0.78(5)(4) 


0.73(3)(7) 
0.66(4)(1) 
0.78(11) 
0.75(12) 


Fermilab-MILC-HPQCD |5] 
APE [22] 
Khodjamirian et al. [51] 
Ball [55] 


Nf = 2 + 1 LQCD 
iV/ = LQCD 
LCSR 
LCSR 




0.68(4)(3) 








B -^TX 


0.27(7)(5) 


0.23(2)(3) 
0.31(5)(4) 
0.26(5)(4) 
0.258(31) 
0.26(4) 
0.26(5) 


Fermilab-MILC [52] 
HPQCD [7] 
APE [22] 
Ball and Zwicky [15] 
Duplancic et al. [T6] 
Wu and Huang 


Nf = 2 + 1 LQCD 
Nf = 2 + 1 LQCD 
Nf = Q LQCD 

LCSR 

LCSR 

LCSR 


B^K 


0.32(6)(6) 


0.331(41) 
0.36(5) 
0.33(8) 


Ball and Zwicky [15j 
Duplancic et al. [T7] 
Wu and Huang [18] 


LCSR 
LCSR 
LCSR 


Bs^K 


0.23(5)(4) 


0.30(4) 


Duplancic et al. [17j 


LCSR 



Table 5: Comparison of the results for /+(0) of the present work with other calculations, 
obtained from lattice QCD (LQCD) simulations or from light-cone sum rules (LCSR) by various 
groups. Where two errors are quoted the first is statistical and the second is the combined 
systematic errors. 



4 Conclusions 

In this article we have presented a lattice QCD calculation of the form factors associated with 
semileptonic decays of heavy mesons. 

We have performed a quenched calculation on a very fine lattice with (3 = 6.6 (a = 0.04 
fm), which allows us to treat the D meson decays in a fully relativistic setup, and to get 
close to the region corresponding to the physical B meson mass. The importance of small 
lattice spacings for heavy-quark simulations has recently become clear in the context of the 
determination of /d^, the decay constant of the Ds meson. In spite of 0(a) improvement, a 
continuum extrapolation linear in seems to be reliable only for lattice spacings below about 



13 



0.07 fm in the quenched approximation [T9l [20] . Depending on the particular improvement 
condition, even a non-monotonous a dependence can appear on coarser lattices. 

In this work we have investigated to which extent the systematic effects caused by lattice 
discretization and long-ranged extrapolations to the physical heavy meson masses may inffuence 
the results of different lattice calculations in which all other sources of systematic errors are 
treated in a similar way. For these reasons, the results of our study can be directly compared 
with those by the APE collaboration in Ref. [22], which reports a very similar calculation on a 
coarser lattice aX (3 = 6.2 (a ^ 0.07 fm) with the same lattice action and currents. Adjusting 
the APE results so that they comply with our procedure for setting the physical value of the 
lattice spacing, we find quite large discrepancies of roughly 2.5cr {D vr) and 2a [D K). If 
we assume that 0{a?) errors are the dominant source of this effect, the difference in the results 
of the two studies suggests an upper limit on the discretization errors of approximately 0.08 or 
slightly above Icr in our numbers for f^^'^{Q) and 0.23 or 3 — 4cr in the APE results. 

It is, however, to be noted that the interpretation of this difference as a mere discretization 
error is somewhat more ambiguous than in the case of the decay constants considered in [19], [20], 
because the momentum transfer adds another parameter that has to be adjusted before the 
comparison can be attempted. The corresponding comparison for B decays can, in addition, 
be undermined by the long-ranged extrapolations in the heavy quark mass and/or q^. These 
results suggest that, for high-precision phenomenological applications, completely reliable rel- 
ativistic lattice calculations of these form factors could require even finer spacings, and that, 
for dynamical simulations at realistic pion masses, this goal might be difficult to achieve in the 
near future. While we believe that the progress in computational power will eventually allow 
one to realize this formidable task, it is fair to say that, for the moment, the less demanding 
approaches which interpolate between the D meson scale and non-relativistic results provide a 
valid alternative. 

Finally, a few words are in order about the general perspective for calculations of the semilep- 
tonic form factors of heavy mesons. Form factors of B decays at small values of the relativistic 
momentum transfer q^ involve a light meson with momentum up to 2.5 GeV in the final state 
and are very difficult to calculate on the lattice, mainly because no lattice effective field theory 
formulation is known for this kinematics that would allow for the consistent separation of the 
large scales of the order of the heavy quark mass, as implemented in the Soft-CoUinear Effective 
Theory. 

Thus one is left with two choices. The first one is to calculate the form factors at moderate 
recoil {m\ — q^ ©(m^AqcD)) using, e.g., the HQET or NRQCD expansion and then to 
extrapolate to large recoil {m?B—q^ ~ 0(m^)) guided by the dispersion relations. The advantage 
of this approach is that the calculations can be performed on relatively coarse and thus not 
very large (in lattice units) lattices. Therefore dynamical fermions may be included, high 
statistical accuracy can be achieved as well as a better control over the chiral extrapolation. 
The disadvantage is that a reliable extrapolation from the q^ > 12 — 15 GeV^ regime accessible 
in this method to = may be subtle. However, this problem may be alleviated by a promising 
new approach, "moving NRQCD" [56], which formulates NRQCD in a reference frame where the 
heavy quark is moving with a velocity v. By giving the B meson significant spatial momentum, 
relatively low can be achieved for lower values of the final state momentum thus avoiding 
large discretisation effects. 

For the particular case of the B ^ -k semileptonic decay the problem of simulating at 
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large recoil can be avoided, at least in principle, since the shape of the form factor f+{q^) 
can be extracted from the experimental data on the partial branching fraction in different 
bins, see, e.g., Ref. [5?]. The normalization can then be fixed by comparison to lattice data 
in the ~ 10 — 20 GeV^ range. This strategy (see Ref. [58] for a detailed discussion) is 
indeed promising and may lead to a considerable improvement in the accuracy of the \Vub\ 
determination from exclusive B decays provided the combined data analysis using the full 
statistics of the BaBar and Belle experiments (~ 4- 10^ BB pairs) becomes available. However, 
for rare decays, such as -B ^ K*li B K*fi~^fi~ etc., which are likely to take the central stage 
at LHCb and super- i? factories, a similar strategy seems to be unfeasible. 

The second choice are simulations with fully relativistic heavy quarks on very fine and large 
(in lattice units) lattices. This procedure is presently bound to the quenched approximation, 
but the benefit is that the extrapolation in the heavy quark mass and potentially also in is 
of much shorter range. It goes without saying that the inclusion of dynamical fermions and 
the approach to the chiral limit are also important problems, but presently it is impossible to 
address all relevant issues within one calculation. 

In our opinion both methods are justified and we have chosen the second option in this paper. 
It turns out that our final results for, e.g., the S — >■ tt decays are consistent with determinations 
based on dynamical simulations and LCSR. This may indicate that the quenching effects are 
rather moderate. From our experience, the main problem that limits the usefulness of this 
approach is the construction of sources for the light hadrons which yield a good overlap with 
states of large momentum. It seems that the presently used sources are not good enough in 
this respect. Improved sources have to be developed if a similar calculation is attempted on a 
larger scale in the future. 
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A Simulation results 

Fig. [6] shows a subset of the form factors fo{q'^) and f+{q'^) that we extracted from our sim- 
ulations, for different combinations of the k values for the heavy and spectator quarks, with 
t^decay product = 0.13472. In Table [6] wc present our results for the fits to the simulation data 
with the Becirevic-Kaidalov parameterization [8] according to Eq. ([6]). This parameterization 
uses the vector meson mass tuh*- Our results for rriH* in lattice units are shown in Table [3 
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■(1-a) 


1//3 


OL 


xVd.o.f. 





13519 


0.13 





13519 


0.775+11 


0.5801?;^ 


0.011- 





77 





13519 


0.129 





13519 


0.7241^9 


0.6041^^ 


0.0911^ 


1 


08 
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13519 


0.121 
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0.66li^ 
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13498 





7421^? 
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0.229l^^ 


1 
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13519 
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13498 





67Atll 


0.714li 


0.337l^i 


1 
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13519 


0.121 
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13498 






36411 


0.8081^6^ 
0.841^0 


0.65liJ 
0.7615*1 
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13519 
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0.670l^^ 


0.338l^^ 


1 
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0.730l^^ 


0.4371^^ 
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0.121 





13472 
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4561^^ 


0.82ll^« 


0.709l^^ 
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13519 


0.115 
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3871^^ 


0.829l^^ 


0.7871^^ 





21 





13498 


0.13 





13519 





78311 


1 VQ 

0.5871^^ 


0.0651^3 





87 





13498 


0.129 





13519 





7401^^ 


0.596l?^ 


0.12l^° 


1 


20 





13498 


0.121 





13519 





5361^1 


0.671}^ 


0.44li^ 





78 
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0.115 





13519 





419l«^ 


0.76l|^ 
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0.13 
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7811^° 


0.613+^f 


0.2001^6^ 


1 


03 





13498 


0.129 





13498 





, 1 Q/1 

' '-'^-29 


0.6271^6^ 


0.256l^^ 


1 


24 





13498 


0.121 





13498 







0. 7231^1 


, 1 no 

0.5571^^ 





56 





13498 


0.115 





13498 





4281^? 


0.765l^^ 


0.6915*1 





36 





13498 


0.13 





13472 





77911 


0.6641^? 


0.325l^^ 


1 


18 





13498 


0.129 





13472 





7261^? 


0.6901^^ 


0.386li 


1 


03 






13498 
13498 


0.121 
0.115 






13472 
13472 






4841^1 
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0.8031^^ 
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0.766li 
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0.13 
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821 3^ 


0.565l^t 


0.053l?^ 





82 





13472 


0.129 





13519 





781+i 


0.560l^^ 


0.094+?° 





97 





13472 


0.121 





13519 





565+^° 


0.661}^ 


0.4211^ 





80 






13472 
13472 


0.115 
0.13 






13519 
13498 






4721,1 
802+^^ 


, 111 

0.721}^ 

1 fin 

0.620l^^ 


0.5911^ 
0.1971^8 






69 
86 





13472 


0.129 





13498 





762+30 


O.625I59 


0.2441q5 


1 


15 





13472 


0.121 





13498 





530+^^ 


0.734+^^ 


0.5491^^ 





68 





13472 


0.115 





13498 





4501^0 


0.764in 


0.66511 





59 






13472 
13472 


0.13 
0.129 






13472 
13472 






8031^6 
7511^1 


0.648l^° 
0.673+i 


0.29011 
0.35311 


1 




01 
80 





13472 


0.121 





13472 





524+^^ 


0.7681^6^ 


0.6161?? 





79 





13472 


0.115 





13472 





42511 


0.823l^? 


0.749l?^ 


1 


07 



Table 6: Results of the fits of the lattice data with Eq. ([6]). Note that at vanishing one has 
/+(0) = /o(0), which is given by cbk ■ (1 — a) (fourth column of this Table). 

The parameters /+(0), a and are then extrapolated to the light-quark mass as described in 
Section [3] and illustrated in Fig. O Finally, the parameters describing the physical form factors 
are obtained through extrapolation of these values to the physical B meson mass (or through 
interpolation to the physical D meson mass), according to the heavy-quark scaling laws, see 
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^hcttvy 


^light 


arriH 


arriH* 


0.13 


0.13519 


0.36811i^ 


0.3949+1^ 


0.13 


0.13498 


0.3762+1^ 


0.40171}^ 


0.13 


0.13472 


0.3867+1?, 


0.41101}^ 


0.129 


0.13519 


0.4060+1? 


0.43001}^ 


0.129 


0.13498 


0.41381}? 


0.43661}^ 


0.129 


0.13472 


0.4240l}g 


0.44581}^ 


0.121 


0.13519 


0.66721}^ 


0.68041^0 


0.121 


0.13498 


0.67431}^ 


0.68681}^ 


0.121 


0.13472 


0.68361}^ 


0.69561}^ 


0.115 


0.13519 


0.83691}! 


0.8460li^ 


0.115 


0.13498 


0.84371}^ 


0.85231}^ 


0.115 


0.13472 


0.85271}^ 


0.86111}^ 



Table 7: Masses of the heavy-light pseudoscalar and vector states, H and H*, respectively, in 
lattice units, for the different {Kheavy, taught) combinations. 

Fig.IZl 
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Figure 6: Sample of form factors directly measured in our simulations. The three panels 
show the results obtained for decay product = 0.13472, and for different values of ^spectator- In 
each plot, the results for /q (denoted by full symbols) and for /+ (empty symbols) are plotted 
against the square of the transferred momentum g^. The results for different values of Hheavy are 
displayed using different symbols: diamonds {i^heavy = 0.13), squares {i^heavy = 0.129), circles 
{i^heavy = 0.121) and triangles {i^heavy = 0.115). 
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Figure 7: Extrapolation and interpolation of the (chirally extrapolated) Becirevic-Kaidalov 
parameters in the heavy quark mass. The plots show the results obtained for the combinations 
/+(0) ■ mj^ (l.h.s. panel), (/3 — l)m// (central panel) and (1 — a)mii (r.h.s. panel), for decays 
of a -B (red squares) or of a D meson (green squares) to a kaon. The results are obtained using 
linear fits in mil (solid lines) according to Eq. ([8]); for comparison, the curves resulting from 
fits to quadratic order in (dotted lines) are also shown. 
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